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We show that the gauge invariant treatment of the Schwinger-Dyson equations of QCD leads to an 
infrared finite gluon propagator, signaHng the dynamical generation of an effective gluon mass, 
and a non-enhanced ghost propagator, in qualitative agreement with recent lattice data. The trun- 
cation scheme employed is based on the synergy between the pinch technique and the background 
field method. One of its most powerful features is that the transversality of the gluon self-energy 
is manifestly preserved, exactly as dictated by the BRST symmetry of the theory. We then ex- 
plain, for the first time in the literature, how to construct non-perturbatively a renormalization 
group invariant quantity out of the conventional gluon propagator. This newly constructed quan- 
tity serves as the natural starting point for defining a non-perturbative effective charge for QCD, 
which constitutes, in all respects, the generalization in a non-Abelian context of the universal 
QED effective charge. This strong effective charge displays asymptotic freedom in the ultravio- 
let, while in the low-energy regime it freezes at a finite value, giving rise to an infrared fixed point 
for QCD. Some possible pitfalls related to the extraction of such an effective charge from infrared 
finite gluon propagators, such as those found on the lattice, are briefly discussed. 
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It has been known for a long time that even though the gluon is massless at the level of the fun- 
damental Lagrangian, and remains massless to all order in perturbation theory, the non-perturbative 
QCD dynamics generate an effective, momentum-dependent mass, without affecting the local 
SU {3)c invariance, which remains intact The existence of this mass is discovered by studying 
the Schwinger-Dyson equations (SDEs) of QCD, in a gauge-invariant framework known as the 
pinch technique (PT) [|l|, To obtain massive solutions it is necessary to include longitudinally- 
coupled massless scalars in the Green's functions, which play a role rather Uke Goldstone excita- 
tions, but do not signal any sort of breakdown of local gauge symmetry, which is preserved. Like 
standard Goldstone bosons these massless scalars do not appear explicitly in the S'-matrix; however, 
they play a crucial role in confinement, furnishing the required long range potential. 

An effective low-energy field theory for describing the gluon mass is the gauged non-linear 
sigma model known as "massive gauge-invariant Yang-Mills" [^, with Lagrangian density 

^MYM = ^Gl^-m^Tr [A^, - g-'U{e)d^,U-\d)]\ (1) 

where = Yi^a^aA"^, the are the SU(3) generators (with TrA^A^ = 2dab), and the N x N 
unitary matrix U{d) = exp [/jAqO"] describes the scalar fields 0„. Note that ^mym is locally 
gauge-invariant under the combined gauge transformation 

A'^ = VA^V-'-g-'[d^V]V-\ U' = U{d')=VU{d), (2) 

for any group matrix V = exp [/iAaft)''(x)] , where (o"{x) are the group parameters. One might think 
that, by employing (^, the fields 0„ can always be transformed to zero, but this is not so if the da 
contain vortices. To use the ^mym in one solves the equations of motion for IJ in terms of the 
gauge potentials and substitutes the result in the equations for the gauge potential. One then finds 
the Goldstone-like massless modes mentioned above. This model admits vortex solutions [^] , with 
a long-range pure gauge term in their potentials, which endows them with a topological quantum 
number corresponding to the center of the gauge group [Za? for SU and is, in turn, responsible 
for quark confinement and gluon screening. Specifically, center vortices of thickness ~ where 
m is the induced mass of the gluon, form a condensate because their entropy (per unit size) is larger 
than their action. This condensation furnishes an area law to the fundamental representation Wilson 
loop, thus confining quarks [[l],|3|]. 

Of course ^mym is not renormalizable, and breaks down in the ultraviolet. This breakdown 
simply reflects the fact that the gluon mass m in (|l]) is assumed to be constant, while the solutions 
of the SDEs reveal that the mass is momentum-dependent, vanishing at large [|l|]. Specifi- 
cally, when studying the SDE for the gluon propagator, A(^^), one looks for infrared finite solu- 
tions, i.e. with A^^(O) > 0. Such solutions may be fitted by "massive" propagators of the form 
^ H^^) = + '«^(^^)> where nP'{q^) is not "hard", but depends non-trivially on the momentum 
transfer q^. When the renormalization-group logarithms are properly taken into account in the SDE 
analysis, one obtains, in addition, the non-perturbative generalization of a{q^), the QCD running 
coupling (effective charge). The presence of m^{q^) in the argument of (x{q^) tames the Landau 
singularity associated with the perturbative j8 function, and the resulting effective charge is asymp- 
totically free in the ultraviolet , "freezing" at a finite value in the infrared. 
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The general picture described above has received spectacular confirmation from a plethora of 
lattice studies, spanning a period of several years: the gluon propagator reaches indeed a finite (non- 
vanishing) value in the deep infrared, as predicted by Cornwall. This rather characteristic behavior 
was already seen in early studies [Q], and has been firmly established recently using large- volume 
lattices, for both SU{2) [g] and SU{3) [|] pure Yang-Mills (no quarks included). 

In this talk we will present recent results from a gauge-invariant study of the coupled gluon- 
ghost system of SDEs [^, yielding an infrared finite gluon propagator and a divergent (but non- 
enhanced) ghost propagator, in qualitative agreement with recent lattice data [Q, |6|] ; this behavior 
has also been confirmed within the Gribov-Zwanziger formalism [||]. As the title suggests, we will 
eventually focus on the issue of the infrared finite QCD effective charge. 

Obtaining an infrared finite result for the gluon self-energy from SDEs, without violating the 
underlying local gauge symmetry, is far from trivial, and hinges crucially on one's ability to devise 
a self-consistent truncation scheme that would select a tractable and, at the same time, physically 
meaningful subset of these equations. Specifically, while in QED the Green's functions satisfy 
naive Ward Identities (Wis), in QCD they satisfy complicated Slavnov-Taylor identities (STIs), 
which involve various composite ghost operators. To see how this complicates the truncation pro- 
cedure of the SDEs, consider the STI of the gluon self-energy 

^^'n^v('?)=o. (3) 

Eq. (^ is without a doubt the most fundamental statement at the level of Green's functions that 
one can obtain from the BRST symmetry; it affirms the transversality of the gluon self-energy and 
is valid both perturbatively to all orders as well as non-perturbatively. The problem is that in the 
SDE governing n^v(^) enter higher order Green's functions, namely the fully-dressed fundamental 
vertices of the theory, which satisfy complicated STIs. Thus, whereas in QED the validity of Eq. (^ 
can be easily seen at the level of the SDE, simply because q^T^{p,p + q) = e ^{p + q) —S^^ (p)] , 
in QCD proving Eq. (||) is very difficult, and requires the conspiracy of all full vertices appearing in 
the SDE. Truncating the SDE naively usually amounts to leaving out some of these vertices, and, 
as a result, Eq. (||) is compromised. Instead, the gauge-invariant truncation scheme [^], based on 
the PT [|T], IP and its correspondence with the background field method (BFM) JlOj ] maintains the 
validity of Eq.(|3|) at every level of approximation. 

The gluon propagator in the covariant gauges has the form iAuviq) = Puv(q)Mq^) + ^ '^^f^ 

where denotes the gauge-fixing parameter, P^v(^) = g^v — QhQv/q^ is the usual transverse pro- 
jector, and, finally, A^^{q^) = q^ + ill{q^), with n^v(^) = P/iv(^)n(^^) the gluon self-energy. The 
full ghost propagator D{p^) and its self-energy L{p'^) are related by iD~^{p^) = p^ — iL{p^). In 
the case of pure (quarkless) QCD, the new SD series |^] for the gluon and ghost propagators reads 
(see also Fig. |l|) 

A [q)V^,[q)- [\+G{q^)Y ' 

iD-\p^) = p2 + iX r'' A^v (^)r^ {p,k)D{p + k) , 

iA^^{q)=X jHflD{k + q)AP''{k)H„y{k,q), (4) 
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Figure 1: The PT-BFM gluon-ghost system. The circles attached to the external gluons denote that, from 
the point of view of Feynman rules, they are treated as background fields. 

where A = g^CA, with Ca the Casimir eigenvalue of the adjoint representation [Ca = for SU {N)], 
and = {2n)^^ J d'^k, with d = 4 — e the dimension of space-time. is the standard (asym- 
metric) gluon-ghost vertex at tree-level, and the fully-dressed one. G{q^) is the g^y component 
of the auxiliary two-point function A^y{q), and the function H^y is defined diagrammatically in 
Fig. [H H(yy is related to the full gluon-ghost vertex by q'^H^y {p,r, q) = —iTy {p,r, q) ; at tree-level, 
^CTv = igav- When evaluating the diagrams (a,) we use the BFM Feynman rules the BFM 
fully dressed three-gluon and gluon-ghost vertices are denoted by F^q.^ and F^. 

Notice a point of paramount importance: due to the Abelian all-order Wis that these two full 
vertices satisfy (for all E,), namely 

q^f^^li=iA-^l{k + q)-iA-^]^{k), q^f ^ = iD-\k + q) - iD-\k), (5) 

one can demonstrate that q^ [{a\ ) + (a2)]^v = and q^ [{a^ ) + (a4)]^v = [[TTl]. Thus, unlike other 
treatments in the literature, within this formalism the transversality of the gluon self-energy, i.e. 
Eq. (||), is preserved at every step, in absolute compliance with the BRST symmetry. 

Next, following standard techniques, we express r^Q,^ and as a function of the gluon and 
ghost self-energy, respectively, in such a way as to automatically satisfy the crucial Wis of Eq. (^); 
failure to satisfy these Wis would invariably compromise the transversality of the answer. The 
Ansatz we will use is 



= + [n„/3 {k + q)- {k)] , F^ = F^ - [L{k + q)- L{k)] ; 

q q 

its essential feature, other than satisfying the aforementioned Wis, is the presence of massless, lon- 
gitudinally coupled pole terms, which are instrumental for obtaining A^^(O) ^ [[l^]. These poles 
are not kinematic but dynamical, corresponding to a composite (bound-state) Goldstone excitation, 
enforcing the local gauge invariance. For the conventional ghost-gluon vertex Fy, appearing in 
the second SDE of (^ we will use its tree-level expression, i.e., Ty ^Ty = —py; this is perfectly 
legitimate, since in this formalism the two ghost vertices, Fy and F^, are different. Finally, for 
Hay we use its tree-level value, H^^] . 

In Fig.^, we show the numerical result for A(^^) renormalized at /X = Mj, = 4.5 GeV, and the 
comparison with the corresponding lattice data of Ref.[^. In the right panel of Fig.^ we present 
the dressing function for the ghost propagator, renormalized at the same point. 
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Figure 2: Left Panel:The numerical solution for the gluon propagator from the SDE (black continuous line) 
compared to the lattice data of Ref.[^]. Right panel:The ghost dressing function p^D{p^) obtained from the 
SDE. In the insert we show the lattice data for the same quantity; notice that there is no "enhancement" [|l3|. 



Let us now turn to the QCD effective charge. There are two main issues: (i) how to define it 
consistently at the level of perturbation theory: specifically, which graphs determine the running, 
and (ii) how to extend the (whatever) definition one reaches in (i) into the non-perturbative regime. 

Point (i) has been addressed exhaustively in the literature [|^] : the upshot is that in the context 
of the PT one may replicate to all-orders in perturbation theory the prototype QED construction of 
an effective charge. To fix the ideas, the PT one-loop gluon self-energy reads 



.2,„ ( q' 



(6) 



where b = \ ICa/^^ti^ is the first coefficient of the QCD j8 -function. Due to the Abelian Wis satis- 
fied by the PT effective Green's functions, the new propagator-like quantity A^^{q^) absorbs all the 
RG-logs, exactly as happens in QED with the photon self-energy. Equivalently, since Zg and Za, the 
renormalization constants of the gauge-coupling and the effective self-energy, respectively, satisfy 
the QED relation Zg = Z^ , the product d{q^) = g-^A{q-^) forms a RG-invariant (jtt -independent) 
quantity [[T]]; for large momenta q^, 

d{q ) = — n— , (7) 
q^ 

where g^{q^) is the RG-invariant effective charge of QCD, 

^^^^^^ " l+bgHn{qyn^) " b\n{qyA^) ' 

Let us now come to point (ii): assuming that one has non-perturbative information about the in- 
frared behavior of the conventional gluon propagator A{q^), how should one extract an effective 
charge, which, perturbatively, will go over to Eq. (||)? To accomplish this, one must use an addi- 
tional field-theoretic ingredient: the conventional A(^^) and the PT-BFM A{q^) are related by the 
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Figure 3: Left panel: The solution of SDE renormalized at jj. ~ — 4.5 GeV (continuous blue curve) and 
ji = Mz = 91 GeV (red line+square curve). Right Panel: The corresponding PT-BFM A{q^) obtained as the 
convolution of A{q^) and the function g'^{n^)/[l + G{q^)]^. 



formal all-order relation [15] 



(9) 



Note that the G{q^) already appears in Eq. (^ and Fig.[I[ With our approximations its SDE reads 



Giq^ 



X 

Jk 



2 + 



Bq^ 



A{k)D{k + q). 



(10) 



First of all, it is easy to verify that at lowest order the G{q ) obtained from Eq. ( |10| ) restores the 
jS function coefficient in front of ultraviolet logarithm. In that limit 1 + G{q^) = 1 + it^t- ln(g^/ }JL^) 



and A- 1(^2) = + f ^ln(^VM^) • Then using Eq. 



4 A?,Tfi 

we recover the A^'(^^) of Eq. (^), 
as we should. Then, non-perturbatively, one substitutes into Eq. (^ the G{q^) and A{q^) obtained 
from solving the system in Eq. (^, to obtain A(^^). This latter quantity is the non-perturbative gen- 
eralization of Eq. (^; for the same reasons explained above, when multiplied by it should form 
an RG-invariant quantity, e.g. the non-perturbative generalization of d{q^). In Fig.^we present the 
combined result of the above steps: d{q^) is obtained from two different sets of solutions of the 
system Eq. (^, one renormalized at ji. = My = 4. 5 GeV and one at /i = Mz = 91 GeV. Ideally the 
two curves of d{q^) should be identical; even though this does not happen, due to the approxima- 
tions employed when solving the system of Eq. (^, the two curves are fairly close, indicating that 
d{q^) is to a very good approximation an RG-invariant quantity, as it should. 

We are now in the position to define the non-perturbative QCD effective charge from the RG- 
invariant quantity d{q^). Of course, given that d{q^) reaches a finite value in the deep infrared, 
it would be completely absurd to define the effective charge by forcing out a factor of \/q^; such 
a procedure would furnish a completely unphysical strong QCD coupling, namely one that would 
vanish in the deep infrared(!) where QCD is supposed to be strongly coupled. 
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Figure 4: The QCD effective charge, 
a[q'^) —'g'^{q'^)/AK, extracted from Fig. 3 by 
factoring out a gluon mass of m(0)=500 MeV. 



The correct thing to do is to factor out a "massive" 
propagator, i.e. write 



(11) 



Of course, as we have emphasized, m^{q^) it- 
self is running, which must also be taken into 
account in a more sophisticated treatment. For 
the purposes of this talk, however, we assume 
that m^{q^) is constant, m^{q^) = m^(0), and use 
for m(0) the value of 500 MeV favored by phe- 
nomenology [16]. The a{q^) obtained is shown 
in Fig.@; as announced, at low energies freezes 
to a finite value, indicating the appearance of an 
infrared fixed point of QCD. 
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